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Summary. Three kinds of magnetic couplings in the quiet solar atmosphere are 
highlighted and discussed, all fundamentally connected to the Lorentz force. First 
the coupling of the convecting and overshooting fluid in the surface layers of the 
Sun with the magnetic field. Here, the plasma motion provides the dominant force, 
which shapes the magnetic field and drives the surface dynamo. Progress in the un- 
derstanding of the horizontal magnetic field is summarized and discussed. Second, 
the coupling between acoustic waves and the magnetic field, in particular the phe- 
nomenon of wave conversion and wave refraction. It is described how measurements 
of wave travel times in the atmosphere can provide information about the topogra- 
phy of the wave conversion zone, i.e., the surface of equal Alfven and sound speed. 
In quiet regions, this surface separates a highly dynamic magnetic field with fast 
moving magnetosonic waves and shocks around and above it from the more slowly 
evolving field of high-beta plasma below it. Third, the magnetic field also couples to 
the radiation field, which leads to radiative fiux channeling and increased anisotropy 
in the radiation field. It is shown how faculae can be understood in terms of this 
effect. The article starts with an introduction to the magnetic field of the quiet Sun 
in the light of new results from the Hinode space observatory and with a brief survey 
of measurements of the turbulent magnetic field with the help of the Hanle effect. 



1 The magnetic field of the quiet Sun 

Over the past three and a half years, the Sun stayed in a minimum state of 
magnetic activity as it has ended cycle 23 and is about to start with cycle 24 
(if not pausing for yet a longer period of time) . In this period of quiescence 
it was possible to observe the Sun with an exceptional instrument, the Solar 



Optical Telescope SOT onboard the Hinode space observatory (Kosugi et al. 



20071. The Japanese Hinode satellite was put in orbit on September 22, 2006. 
It is not so much the spatial resolution of 0.3" that makes this instrument 
exceptional for quiet- Sun observing, but rather the absence of seeing in com- 
bination with high pointing accuracy. This allows for unprecedented "deep" 
(long-exposure) polarimetry with correspondingly high polarimetric sensitiv- 
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ity at a spatial resolution of 0.3". A similar polarimetric accuracy at this high 
spatial resolution has not been achieved from the ground in the past. 

Immediately evident in the total or circular polarization maps over a large 



field of view from Hinode (Lites et al. 2008) is the magnetic network, which 
persisted existing during the most quiet states of the Sun. The magnetic 
network consists of an accumulation of magnetic fields in the borders between 
supergranular cells. The origin of these magnetic fields remains an enigma. 
Is it the result of advection by the supergranular flow to which it is often 
described, or, vice versa, is the supergranular flow rather a consequence of the 
existence of the magnetic network (R. Stein, private communication)? Is the 
network fleld generated locally, near the surface, or is it an integral part of 
the globally acting dynamo ( Stein et "ar]|2003 ), or is it just the decay product 
of sunspots and/or ephemeral active regions? 

Also omnipresent in this quiescent state of the Sun are small-scale mag- 
netic field concentrations, visible as delicate, bright objects within and at 
vortices of intergranular lanes. The structure made up of ensembles of bright 



elements is known as the filigree (Dunn & Zirker 1973). Mehltretter (19741, 



while observing in the visible continuum, referred to them as facular points 
because they are the footpoints of magnetic field concentrations that appear 
as faculae near the solar limb. In more recent times, these objects were mostly 



observed in the G band (a technique originally introduced by MuUer 1985) 



because the molecular band-head of CH that constitutes the G band, acts as 



a leverage for the intensity contrast 


(Rutten|1999 


Rutten et al.|2001 
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1. Being located 



in the blue part of the visible spectrum, this choice also helps improving the 
diffraction-limited spatial resolution and the contrast in the continuum. Re- 
cent observational investigations of the dynamics, morphology and properties 
of small-scale magnetic field concentrations of the quiet Sun includ e [Berger 
et all (|2004P; ILanghans et al.l (|2004|); ILites & Socas-Navarrol (|2004|) ; ISanchez 



Almeida et al. (2004); Socas-Navarro & Lites (2004); Wiehr et al. (20041; 



Rouppe van der Voort et al. (2005); Dommguez Cerdena et al.| (| 2006a|b l; 
Berger et all (|2007|); iBovelet & Wiehrl (|2007|) ; iCenteno et al.l (12007) ; Ishikawa 



et al. (2007); Langangen et al. (2007); Martmez Gonzalez et al.|(,2007|);|Orozco 
Suarez et al.| (|2007 P; |Rezaei et al.| (|2007a|b|); [Tritschler et al.| (|2007 ~|Bello 



Gonzalez fc Kneer (2008); Bcllo Gonzalez et al. (2008); de Wijn et al 



(2008); 



Qrozco Suarez et al.| ( |2008| ); |Bello Gonzalez et al.| ( |2009[ ) and references therein. 
It would be interesting to learn to which degree the abundance of small-scale 
magnetic field concentrations persists in a grand minimum, to find out more 
about their origin. 

With the help of the spectropolarimeter of the Solar Optical Telescope 
(SOT) onboard Hinode it became for the fist time possible to reliably deter- 
mine the transversal (with respect to the line-of-sight) magnetic field compo- 
nent of the quiet Sun. These measurements indicate that, seen with a spa- 
tial resolution of 0.3", the quiet internetwork regions harbor a photospheric 
magnetic field whose mean field strength of its horizontal component con- 
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siderably surpasses that of the vertical component (Lites et al. 2007 2008 



Orozco Suarez et "ar]|2007 1 . According to these papers, the vertical fields are 



concentrated in the intergranular lanes, whereas the horizontal fields occur 
most commonly at the edges of the bright granules, aside from the vertical 
fields. Lites et al. ( 2008 1 determine for the horizontal field component a mean 
apparent field strength (averaged over a large field of view including network 
and internetwork regions) of 55 G, while the corresponding mean absolute ver- 
tical field strength was only 11 G. Harvey et al. (20071 find from recordings 
with GONG and SOLIS at moderate angular resolution a "seething magnetic 
field" with a line-of-sight component increasing from disk center to limb as 
expected for a nearly horizontal field orientation. It is reasonable to assume 



that the horizontal fields of Lites et al. (2008) and those of Harvey et al. 



(2007 



are different manifestations of the same magnetic field. Ishikawa et al. 



( 2008 1 detected transient horizontal magnetic fields in plage regions as well. 
Previously, Lites et al. (19961) andlMeunier et al. (1998) reported observations 
of weak and strong horizontal fields in quiet Sun regions. 

The anisotropy of the quiet-Sun magnetic field as revealed by the Hinode 
measurements is not necessarily in contradiction with the observed depolar- 
ization of scattered light through the Hanle effect. For the quantitative in- 
terpretation of the Hanle effect it was customarily assumed in the past that 
the quiet Sun magnetic field was "turbulent" in the sense that it was isotrop- 
ically distributed on the scales relevant for the analysis. Theories have now 
been developed to include distribution functions for the field strength and 



angular directions (see, e.g., Carroll & Kopf 2007 Sampoorna et al. 2008 



Sampoorna 2009 ) for a better description of the turbulent magnetic field in 



radiation transfer and first steps are taken to use such formulations for the 



computation of the Hanle effect in scattering media (Frisch 2006 Nagendra 

eFaTlpoogl 



The isotropy assumption may still be valid on scales smaller than 0.3". In 
fact, simulations suggest constant angular distribution within ±50° from the 
horizontal direction on a scale of 0.05". On the other hand, should the field 
be strongly anisotropic on all scales, it would still produce Hanle depolariza- 
tion but its interpretation would be less straightforward. From a theoretical 
point of view, the anisotropy of the magnetic field comes not as a surprise, 
and homogeneous turbulence of the magnetic field seems unlikely since the 
convective flow is far from homogeneously turbulent given the scale size of 
1000 km of granules vs. the pressure scale height of 100 km. 



Stenflo ( 1982 ) roughly estimated the strength of the weak turbulent mag- 



netic field based on the Hanle technique between 10 and 100 G. This range 
was narrowed to 4-40 G by Stenfio et al. (1998). More precise estimates by 



Faurobert-SchoU (1993) and Faurobert-SchoU et al. (1995) yielded values in 



the range of 30 to 60 G in the deep photosphere and 10 to 20 G in the middle 
and upper photosphere, where the latter values were increased to 20 to 30 G by 



Faurobert et al. (2001 ). The values reported by Stenfio et al. and by Faurobert 



et al. are not without controversy, however. Trujillo Bueno et al. (2004) find 
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from three-dimensional radiative transfer modelling of scattering polarization 
in atomic and molecular lines an ubiquitous tangled magnetic field with an 
average strength of 130 G around 300 km height in the photosphere, which is 
much stronger in the intergranular than in the granular regions. They estimate 
that the energy density of this field would amount to 20% of the kinetic energy 
density of the convective motion at a height of 200 km. If this high value is 



correct it also indicates that the Zeeman measurements with Hinode (Lites 



et al. 2007 2008 Orozco Suarez et al. 20071 have not captured quite all of 



the existing quiet-Sun fields, presumably because of polarimetric cancellation 
which Zeeman measurements are subject to in contrast to Hanle measure- 
ments. Adopting the idealized model of a single-valued microturbulent field, 



Trujillo Bueno et al. ( 2004 1 obtained a mean field strength that varies between 



50 to 70 G in the height range from 400 km to 200 km, respectively, and only 
when taking an exponential probability density function for the field-strength 
distribution into account do they obtain the higher value of 130 G. Relatively 



(20051. 



high values are also reported by Bommier et al 

More in line with Faurobert et al. are Shapiro et al.| (|2007|, who ob- 
tain from differential Hanle measurements with the CN violet system a field 
strength in the range from 10 to 30 G in the upper solar photosphere, while 
the analyses of the observed scattering polarization in C2 lines by [Faurobert] 
k. Arnaud (2003) and Berdyugina & Fluri ( |200 4| imply a field strength of 



about 10 G. However measurements of the scattering polarization in molecu- 
lar lines may be quite sensitive to the thermal structure in the atmosphere. 
In fact Trujillo Bueno et al. (2004) obtain from measurements with C2 a field 
strength of the order of 10 G too, but they also show that these measurements 
sample the atmosphere mainly above granules only, where, correspondingly, 
the turbulent field must be much weaker than in the downflows of the inter- 
granular space (see the review by Trujillo Bueno et al. 2006 for a detailed 
presentation of their Hanle-effect measurements). 

In any case, it seems that "deep" polarimetric measurements with Hinode 
have discovered a large part of the hitherto "hidden" magnetic field that was 
known to us only through Hanle measurements. It made this field accessible 
to Zeeman analysis and therefore to a more reliable determination of its angu- 
lar distribution, at least down to 0.3" spatial resolution. In the next chapter, 
we review results from recent simulation that aim at explaining the predom- 
inance of the mean horizontal over the mean vertical field in the quiet-Sun 
photosphere. 



2 Coupling of convection with magnetic fields 

It was mentioned in the previous chapter that a few observational studies 
prior to Hinode already hinted at a frequent occurrence of horizontally ori- 
ented magnetic fields in the quiet Sun. Likewise, the horizontal fields did 



not come unannounced to theoretical solar physics. Grossmann-Doerth et al. 
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( 1998 1 noted "we find in all simulations also strong horizontal fields above 



convective upflows", and Schaffenberger et al. (2005 20061 found frequent 



horizontal fields in their three-dimensional simulations, which they describe 
as "small-scale canopies". Also the three-dimensional simulations of |Abbett] 
(20071 display "horizontally directed ribbons of magnetic flux that perme- 
ate the model chromosphere" , not unlike the figures shown by [Schaffenberger 
et al. (2006). However, these reports did not receive wide attention because 



actual measurement of the weak transversal component was not possible or 
unreliable prior to the advent of Hinode. 
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Fig. 1. Left: Mean absolute horizontal magnetic field components, {|-Bi|) (— ■ — • 

— ■ — ), and {\By\) ( ), and absolute vertical field component, (|-Bver|) ( ) 

as a function of optical depth logrgsonm of the dynamo run of [Schiissler fc Vogler| 
(20081. The averaging refers to surfaces of constant reaonm- The solid curve is the 

1/2 



and (\Bv 



(- 



rms {Bl + By}"'\ Right 
as a function of height z from the simulation run h20 (heavy) and run vlO (thin) of 
(2008]). vlO and h20 substantially differ in their initial and boundary 
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conditions for the magnetic field. Note the different physical meanings of the abscissa 
and the different units and scales in the ordinates of the two plots. 



More recently and after the discovery of the horizontal field with Hinode, 
two theoretical works ( Schiissler & Vogler 2008 and Steiner et al.|2008 1 specif- 
ically aimed at finding out more about its nature and origin. Both papers 
present results of three-dimensional magnetohydrodynamic numerical sim- 
ulations of the internetwork magnetic field with regard to the intrinsically 
produced horizontal magnetic field. In the following I briefly summarize and 
compare part of their results. 

The two simulation runs presented by 



Steiner et al. (2008) and the "local 



dynamo run" of |Schussler fc V5gler (2008) differ substantially in their initial 
and boundary conditions for the magnetic field. Yet, they all show a clear 
dominance of the horizontal field in parts or the full height range where the 
spectral lines used for the Hinode observations are formed. Thus, the intrinsic 
production of a predominantly horizontal magnetic field in the photosphere of 
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three-dimensional magnetohydrodynamic simulations is a rather robust result. 
Figure [T] shows the horizontal and the vertical magnetic field strengths as a 
function of height in the atmosphere of the three simulations. Left and right 
boundaries of the left panel correspond to approximately z — 400 km and 
z = —1000 km of the right panel, respectively. Note that the scale of the 
ordinate is logarithmic and in gauss in the left panel but linear and in mT in 
the right panel. Also account for the non-linear relation between the abscissa 
of the two panels. Interestingly, both simulation runs of Steiner et al. ( 2008 ) 
show a local maximum of the horizontal field component near 500 km height 
and this is also the case for a local dynamo run when the top (open) boundary 
is located at z = 650 km (M. Schiissler private communication) . 

How do these results compare with Hinode? For a fair comparison it is 
indispensable to synthesize the Stokes profiles of the 630 nm Fe I spectral line 
pair from the simulations and subsequently derive whatever parameters were 
derived from the actual observations. The analysis of the synthetic data must 
proceed in the very same manner as done with the observed profiles. Applying 
the appropriate point spread function ( Wedemeyer-Bohm 20081 to the syn- 
thetic profiles and subjecting them to the same procedure for conversion to 
apparent fiux densities as done by Lites et al. (2008) for the observed profiles. 



Steiner et al. (20081 obtain spatial and temporal averages for the transver- 
sal and longitudinal apparent magnetic flux densities, \B^^J, and \B, 



'appi 



'appi 



of 



respectively 21.5 G and 5.0 G for run h20 and 10.4 G and 6.6 G f or run vlO. 
Thus, the ratio r = {\Bf 

obtain from Hinode SP data 



(|2008') 



appl)/(|SaVl> = 4.3 for h20 and 1.6 for vlO. 

app I 



Lites et al. 



11 G re- 



.,pp„ - 55Gand (iSi, 
suiting in r = 5.0. Run vlO was judged to rather reflect network fields because 
of its preference to produce vertically directed, unipolar magnetic fields, en- 
forced by its initial and boundary conditions. For the internetwork field, h20 
is more appropriate. Correspondingly, the r-value of h20 better agrees with 
the measurements of Lites et al. (2008), which measures mainly internetwork 
magnetic fields. It should be cautioned that r is quite dependent on spatial 
resolution in the sense that lower resolution overestimates this value. The rea- 
son for this behavior is that horizontal fields have a more patchy, smoother, 
and less intermittent character than the vertical fields and are therefore less 
subject to polarimetric cancellati on |^ For a com parison between synthetic and 



observed center-to-limb data see Steiner et al. (20091 



What kind of physical process produces the horizontal fields? [Schiissler fc| 



Vogler ( 2008 ) and Steiner et al. ( 2008 1 offer two different but not necessarily 



exclusive explanations. Rather they emphasize two different aspects of the 
coupling of convection with magnetic fields that is at the origin of the hori- 



^ Remember that polarimetric cancellation also occurs for independent field com- 
ponents of the horizontal field when they are perpendicular to each other within a 
single pixel area. Polarimetric cancellation occurs for horizontal fields in the same 
way as for vertical fields - transversal perpendicular fields lead to polarimetric 
cancellation just as antiparallel longitudinal fields do. 
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Fig. 2. Flux expulsion in a close-up from a MHD simulation by Schaffenberger et al 
| |2005[ ): Logarithmic magnetic field strength in a vertical cross-section (top) and in 
three horizontal cross-sections (bottom) at heights of km, 250 km, and 500 km. 
The emergent intensity is displayed in the rightmost panel. The arrows represent 
the velocity field in the shown projection planes. The white curve in the upper panel 



marks the height of optical depth unity. From Wedemeyer-Bohm et al. ( 2008 \ 



zontal fields. Steiner et al. (20081 emphasize the aspect of the flux expulsion 



process ( Weiss||1966 Galloway fc Weiss||1981 1, which describes the expulsion 
of magnetic field from the interior of an eddy fiow like that of granules. Thus, 
the fact that the magnetic field tends to be located in the intergranular space 
and not within granules is considered a consequence of the flux expulsion 
process. However it should be noted that the granular flow is not bounded 
alone by intergranular lanes but also by the overlaying photosphere, which 
efficiently damps overshooting flow owing to its superadiabatic stratification. 
Hence, magnetic field tends not only to be expelled in the lateral direction to 
the intergranular lanes but also in the vertical direction, where it accumulates 
in the upper photosphere and lower chromosphere. In fact, vertical sections 
through the computational domain such as Fig. 1 of Schiissler &: Vogler ( 2008 ) 
and Fig. 3 of Schaffenberger et al.^(2005) show magnetic voids where the gran- 
ular flow is most vigorous as a consequence of the flux expulsion process. The 
voids are arched by horizontally directed magnetic field. This can be nicely 
seen in the close-up shown in Fig. [2] 
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Fig. 3. Dynamo action amplifies and maintains tlie magnetic field in the convec- 
tively unstable layer below Tc = 1. In the layers above, the field is mainly determined 
by its distribution at the Tc = 1 surface. This configuration leads to a steep decline 
of the absolute vertical flux with height as can be seen by counting the loop foot- 
points at each indicated level because small loops are more abundant than large 
loops. On the other hand, this configuration leads to a less steep decline of the mean 
horizontal field strength. 



Schiissler & Vogler ( 2008 1 emphasize the aspect from the local dynamo 



that operates in the convectively unstable layers beneath the surface of 
continuum optical depth unity. Near the surface, weak magnetic field gets 
quickly stretched and thus amplified by the convective flow, in particular also 
by the small-scale turbulent flow of intergranular downflowsj^ On the other 
hand, in the convectively stable photosphere above, the flows become weaker 
and field amplification rapidly drops with height. Thus, the magnetic field in 
the photosphere and its decay with height is mainly determined by the field 
distribution at the surface Tc = 1, in particularly by its energy spectrum as 
a function of horizontal wave number, which in turn is determined by the 
turbulent dynamo beneath this surface. This results in a steep decline of the 
absolute vertical magnetic fiux with height as can be seen from Fig. [3j While 
many loops of small scales (where the energy spectrum is maximal) contribute 
to the vertical flux in the deep photosphere, fewer loops of large scales (where 
the energy is less) add to it in the higher layers. On the other hand, this 
configuration leads to a less steep decline of the horizontal field and hence, 
the mean horizontal field strength starts to dominate the mean absolute ver- 
tical field strength as a function of height. In this picture the dominance of 
the horizontal field component is a natural outcome of the dynamo-generated 
field. 



There is nothing mystic about this amplification, which is a natural consequence 
of the field being tied to the plasma in (quasi) ideal MHD, which does work 
against the Lorentz force on the expense of kinetic energy. However, |Vogler ^ 



Schiissler ( 2007 1 were able to demonstrate that a local dynamo operates in these 



layers, which means that a magnetic field of constant mean energy density is 
maintained without the need of continuous supply of a weak (seed) field. It even 
survives when downfiows continuously pump magnetic field out of the simulation 
box. 
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h = 0.1-L= 100 km 




a.L = a. 1000 km 




Fig. 4. Two intergranular magnetic flux concentrations of opposite polarity are 
connected by a small-scale loop (canopy) in the photosphere, spanning a distance 
of granular scale. Even though this magnetic field configuration is divergence free, 
it follows that the measured mean absolute horizontal field is much stronger than 
the mean absolute vertical field: {|Bh|)/(|_Bv|) ~ 5a. 

It was argued in the course of this conference by J.O. Stenflo that a pre- 
dominance of the horizontal field over the vertical one was in contradiction 
with the solenoidality condition for the magnetic field. Leaving aside that the 
simulations strictly maintain solenoidality and still show a predominance of 
the horizontal over the vertical component. Fig. [4] provides another counter 
example to this conjecture. Assume that within an area A of granular size 
L there are two vertical flux concentrations of, say i3v500G, that occupy an 
area of fvA, where /v ~ 0.02. The two flux concentrations of opposite polarity 
are connected by a photospheric, solenoidal arch of thickness h ~ O.IL as it 
occurs in simulations (viz. the "small-scale canopies" of [Schaffenberger et al.| 
[2005 J . Then flux conservation demands that in a cross section of the loop (as 
the one indicated by the dashed line) <P — B\^Lh — BvfvA/2, where is the 
vertical field strength at Tc w 1 and the strength of the horizontal field 
of the arch. It follows that B^ — ^taf^B^. With a w 1 we obtain B^ w 50 G. 
When the horizontal field of the arch fills about f}^ — 0.8 of the area A, we 
obtain a mean horizontal field of (|-Bh|) = ^h/zi = B^fh/v^a = 40G while 
(|i?v|) = B^fv = 10 G. The ratio (|i3h|)/(|-Bv|) = 5a/h can be made arbitrarily 
large by increasing a, i.e., by stretching the arch. 

Confusion may arise because of misinterpreting transversal Zeeman mea- 
surements as being a measure of magnetic fiux. It is true that for longitudinal 
Zeeman measurements, the measured mean fiux density is directly propor- 
tional to the magnetic fiux (disregarding polarimetric cancellation effects for 
now). Such a generalization is not valid for the transversal Zeeman effect, 
where the mean flux density is an average value over an area that is not 
perpendicular but parallel to the measured field component. In this case the 
apparent fiux density is a spatially averaged fiux density, not more. Assuming 
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that the horizontal field has the vertical extent of the atmospheric scale height, 
h, with a filling factor of 1, then the horizontal magnetic flux is — {Bh)hL, 
where L is the scale of the field of view. The vertical flux is (p^ ~ {Bv)L'^ 
and consequently, using flux conservation (Bi^) / {B^) oc L/h. On granular 
scales L/h w 1000 km/100 km ~ 10, which again suggests that the predomi- 
nance of the apparent horizontal flux density over the vertical one is a direct 
consequence of the anisotropic nature of convective turbulence in the highly 
stratified atmosphere of the Sun, where the granular scale is ten times larger 
than the pressure scale height. However, remember that this calculation does 
not apply to subgranular scales at which the field may still be isotropic. 

3 Coupling of waves with magnetic fields 

Acoustic waves are generated and emitted by the convectively unstable layers 
beneath the solar surface (rsoonm = !)■ They couple to the magnetic field 
when they enter the atmospheric layers above. Such coupling can be seen to 
take place in Fig. [5] which shows a time instant of a two-dimensional simu- 
lation. Starting from a state where a magnetic flux concentration has formed 
at approximately x — 4150km (left panel), a plane-parallel, monochromatic 
acoustic wave is introduced at the bottom of the computational domain (right 
panel) , which propagates within a time span of about 200 s across the non- 
stationary atmosphere - from the convection-zone, z G [—1200, 0] km, through 
the photosphere, z € [0, 500] km, into the magnetically dominated chromo- 
sphere, z G [500, 1600] km, where it gets partially refracted and reflected by 
interaction with the magnetic field. 

The perturbation of the wave front in the convection zone in Fig. |5] is 
not due to the presence of a magnetic field, but rather to the vigurous inter- 
granular downflows and associated temperature deficit. However, as soon as 
the wave front enters the magnetically dominated atmosphere where /? < 1, 
coupling with the magnetic field kicks in and part of the wave gets primarily 
magnetically driven. (/3 is the ratio of magnetic to thermal pressure.) The 
major effects of this interaction are that (i) the wave front speeds up since the 
Alfven velocity becomes the characteristic speed, which sharply increases with 
height when magnetic pressure drops less quickly than the gas pressure, and 
(ii) the waves refract because of the inhomogeneous magnetic field, defining 
an inhomogeneous refractive index for the magneto-acoustic wave. 

When measuring the wave travel time between two fixed geometrical height 
levels in the atmosphere (representing the formation height of two spectral 
lines) we observe a decrease at locations of strong magnetic field concentra- 
tions because of (i). Moreover, due to (ii) the wave behaves like evanescent 
because of the strong refraction that effectively leads to a reflection of the 
wave. Therefore, the wave travel time betrays the presence of the magnetic 
field concentration and it can be used to map the topography of the magnetic 
field in the solar atmosphere. In fact, this effect was employed by [Finsterle 
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Fig. 5. Left: Logarithmic absolute magnetic flux density in a two-dimensional sim- 
ulation domain. Magnetic flux concentrations form in the downdrafts of convection. 
A particularly strong one has formed near x = 4150 km. Right: A plane-parallel wave 
of frequency 20mHz travels through the convecting plasma into the magnetically 
structured photosphere and further into the low-/3 (magnetically dominated) chro- 
mosphere. The panel shows the difference in absolute velocity between the perturbed 
and the unperturbed solution 168 s after launching the wave. The magnetic field at 
launch time corresponds to that of the panel to the left. Optical depth Tsoonm = 1 
is close to z — 0. The velocity scaling is logarithmic with dimension [cm/s]. At the 
location of the magnetic flux concentration the initially fast (acoustic) wave has 
converted character to fast magnetic and it underwent refraction to such a degree 
that the wave front extending from (x, z) — (2400, 1500) to (3400, 500) has already 
completely turned around and travels back into the atmosphere again. A similar 
fanning out of the wave front starts to occur around x = 1100 km. Adapted from 



Steiner et al. (20071 courtesy of Ch. Nutto. 



et al. (2004) in order to obtain the three-dimensional topography of the "mag- 
netic canopy" in and around active regions by determining the travel time of 
high-frequency acoustic waves in the solar chromosphere. Steiner et al. (20071, 
from where Fig. [5] is derived from, demonstrated with the help of numerical 
experiments that wave travel-times can indeed serve this purpose. 

The theory and theoretical aspects of magneto-acoustic waves in a grav- 
itationally stratified atmosphere, sometimes called magneto-acoustic-gravity 
(MAG) waves or magneto-atmospheric waves, have received much attention in 



recent years. Pioneering works include those of Thomas ( 1982 ) and Zhugzhda 
& Dzhalilov (1982 1984a|b|c I. Shibata (19831 carried out initial numerical 



computer experiments with magneto-atmospheric waves. In more recent times. 



Rosenthal et al. ( 2002 1 and Bogdan et al. ( 2003 1 published two comprehensive 



papers on the subject. These works include several numerical experiments with 
non-uniform magnetic field equilibria in a two-dimensional, stratified atmo- 
sphere. They recognized and highlighted the role of refraction of fast magnetic 
waves and the role of the surface of equal Alfven and sound speed as a wave 
conversion zone. Aiming at applications in local helio-seismology, Cally (2005) 
derives gravito-magneto-acoustic dispersion relations and then uses these to 
examine how acoustic rays entering regions of strong field split into fast and 
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slow components and the subsequent fates of each. Cally ( 2007 1 presents the 
theory in a particularly instructive manner. Results from numerical simula- 
tions of MAG-wave propagation in three-dimensional space are presented by 



Cally & Goossens (2008) and Moradi et al. (2009) 



Khomenko & CoUados ( 2006 1 carried out numerical simulations of magne- 



to-acoustic wave propagation in sunspots and found that the fast (magnetic) 
mode in the region where Cs < va does not reach the chromosphere but 
reflects back to the photosphere due to wave refraction, caused primarily by 
the vertical and horizontal gradients of the Alfven speed. For small-scale flux- 
tubes, Khomenko et al. ( 2008b I find that deep horizontal motions of the flux 
tube initially generate a slow (magnetic) mode and a surface mode that are 
efficiently transformed into a slow (acoustic) mode when the magnetic field 
starts to dominate. This slow mode propagates vertically along the magnetic 
field remaining always within the flux tube, where it steepens to a shock. 
Only a small part of the driver energy is returned to the photosphere by 
the fast magneto-acoustic mode. Khomenko et al. (2008a) demonstrate that 



photospheric five-minute oscillations can leak into the chromosphere inside 
small-scale vertical magnetic fiux tubes as a consequence of radiative damping, 
which leads to a significant reduction of the cutoff frequency and they provide 
observational evidences of this effect. This effect is not to be confounded with 
the "ramp effect" (Cally 2007), which lowers the cutoff frequency when the 



flux tube is inclined with respect to the gravitational acceleration (Suematsu 
1990 Jefferies et al.||2006 l. Both these works of Khomenko et al. suggest that 



vertical magnetic field concentrations play an essential role in coupling the 
dynamics of the photosphere to the chromosphere through efficient channeling 
and conversion of magneto-acoustic waves. 

Figure [6] demonstrates the complexity of magneto-acoustic wave propa- 
gation in a magnetically structured, stratifled atmosphere. A magnetic flux 
sheet (two-dimensional) of a strength of 1600 G at its base (where /3 < 1) is 
shifted to the right in the transverse direction with a single impulse of 12 s 
duration and a maximal velocity of 0.75 km/s. As a consequence of this sud- 
den movement, four types of waves emanate from the base of the flux sheet: 
(i) a fast (magnetic) wave in the low-/3 regime of the flux-sheet interior, visi- 
ble in the middle panel, (ii) a slow (acoustic) wave that propagates along the 
magnetic fleld within the flux sheet, visible in the left panel, (iii) a fast acous- 
tic wave that propagates spherically into the ambient medium, visible in the 
right panel, and (iv) a slow magnetic wave that propagates in the boundary 
layer of the flux sheet. The predominantly acoustic waves show an antisym- 
metric pattern with respect to the sheet axis because of the movement to the 
right, which causes a compression at the leading edge and a decompression 
at the trailing edge of the flux sheet. However, as the slow wave propagates 
it becomes asymmetric (not visible in Fig.|6|. On the left side the wave crest 
quickly steepens into a shock because of being preceded by a wave trough. On 
the right side the inverse sequence causes the wave to spread, which impedes 
at first the development of a shock. 
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Fig. 6. Snapshot of magneto-acoustic waves, generated by an initial impulsive 
movement of the equilibrium magnetic flux concentration (black field lines) to the 
right. Left: Longitudinal velocity of the slow (predominantly acoustic) wave that 
propagates along the magnetic field. Middle: Transversal velocity due to the fast 
(predominantly magnetic) wave within the flux sheet. It undergoes differential re- 
fraction. Right: Temperature perturbation due to these two waves and due to the 
fast acoustic wave in the ambient medium. The white contours indicates f3 — 1. 
Velocities are only shown where \B\ > 50 G. Courtesy G. Vigeesh. 



The shape of the fast (magnetic) wave becomes crescent (middle panel) 
because of the non-uniform Alfven speed. As the tips of the crescents move 
essentially sideways, they enter the low-/3 periphery of the flux sheet where 
they convert to predominantly acoustic (fast) and therefore become visible in 
ST (right panel) as the wing like feature that extends in the vertical direction. 
From a rough estimate of the acoustic energy flux generated by such impulsive 
transverse motions, Vigeesh et al. ( 2009 1 conclude that this flux would hardly 
balance the chromospheric energy requirements in the network. Previously, 
similar numerical experiments with impulsive and periodic driving have been 
carried out by Hasan et al. (20051 and Hasan & van Ballegooijen (20081. 

Until here we have considered the coupling of waves with intense, small- 
scale magnetic flux concentrations like they occur in network and plage re- 
gions. The numerical experiments mentioned above, all make use of highly 
idealized wave drivers: typically a monochromatic transversal or longitudinal 
periodic motion, or a single impulsive motion is imposed on an initially static 
equilibrium configuration. Certainly for the internetwork, such quasi static 
states are unrealistic. Furthermore, magnetic flux concentrations in the inter- 
granular lanes of the internetwork attain typically hectogauss not kilogauss 
field strength and they rather connect with the omnipresent horizontal field 
than vertically extending into the chromosphere. In this case, the strength of 
the internetwork field can be expected to exponentially decrease with height 
like the gas pressure and the density do, so that the magnetic field would 
not become dominant in the upper layers of the photosphere and the chro- 
mosphere and consequently no substantial coupling between waves and the 
magnetic field would occur in these layers. On the other hand if there is a pre- 
dominance of one magnetic polarity, part of the magnetic flux can be expected 
to connect to the outer solar atmosphere or to a region of opposite polarity 
further away, in which case the field strength would decrease less steeply lead- 
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ing to /3 < 1 in the chromosphere (but not yet in the photosphere). How do 
waves couple to the magnetic field under these circumstances? 




Fig. 7. Three horizontal cross sections through a a simulation domain with 9.6 Mm 
side length displaying the temperature. Left: Section at z — —200 km showing gran- 
ules (hot) and intergranular lanes (cool). Middle: Section at z — +200 km showing 
inverse granulation. Right: Section at 2 = +1000 km showing the "fluctosphere" con- 
sisting of a rapidly changing network pattern of hot material compressed in traveling 
shock waves. Pockets of cool, expanded material reside in between the hot plasma. 



First to the waves. Figure [7] shows the temperature in three horizontal 
cross sections through the computational domain of 9.6 Mm side length. At 
200 km below the rsoonm = 1 surface we see the cool intergranular lanes 
and hot granules (left panel), at +200 km the hot intergranular lanes and 
cool granules of the inverse granulation (seen a bit higher up in spectroscopic 
quantities). Both these patterns evolve on roughly the granular time scale. In 
the cross section at +1000 km we see a totally different pattern that evolves 
on a much shorter time scale and is due to shock waves that have formed at 
this height range and travel in all directions, forming a network of hot ma- 
terial. This shock- wave pattern that emerges from acoustic waves, which are 
generated by the convective granular motion at the base of the atmosphere 
was first shown to exist in three-dimensional simulations without magnetic 
field by Wedemeyer et al. ( 2004 1 . It leads to large fluctuations in the tenuous 



atmosphere above the classical temperature minimum, to a veritable "flucto- 
sphere" ( Wedemeyer-Bohm et al.|2008 l, earlier termed the "clapotisphere" by 



Rutten & Uitenbroek (19911 and Rutten (19951 with regard to peak fluctua- 



tions caused by the interference of (shock-)waves. 

In combination with magnetic fields, these disturbances give rise to a rich 
variety of magneto-acoustic wave phenomena. As detailed above, the gas pres- 
sure in the gravitationally stratified atmosphere may drop more quickly with 
height than the magnetic energy density does, giving rise to a height range 
where sound speed and Alfven speed are of similar magnitude. Within this 
region, which forms a corrugated surface excursive over a wide height range 
in the three-dimensional atmosphere, propagating wave modes change nature 
from acoustic to magnetic and from slow to fast and vice versa. Above this 



Magnetic Coupling in the Quiet Solar Atmosphere 15 

surface there is a predominant tendency for magnetic modes to get refracted 
and reflected due to the dispersive nature of the inhomogeneous magnetic 
field. 



Schaffenberger et al. ( 2005 1 have simulated this case with a field of a pre- 
dominant polarity of constant mean net vertical flux density 10 G. Their three- 
dimensional simulation domain encompasses a height range from —1500 km 
to -t- 1500 km (where zero corresponds to rsoonm = !)• Immediately apparent 
from a movie that shows the field strength (http : //www . kis . uni-f reiburg . | 
|d.e/~steiner/vsec .mov) is that the surface of /3 = 1 (where Cg « va) sep- 
arates a region of highly dynamic magnetic fields with fast moving magne- 
tosonic waves and shocks around and above it from the more slowly evolving 
field of high-beta plasma below it. This surface is located at roughly 1000km 
in this case. It is corrugated and its local height strongly varies in time over 
a range of about 1000 km. 

The magnetic field in the chromosphere of this simulation continuously 
rearranges itself on a time scale of less than 1 min, much shorter than in 
the photosphere or in the convection-zone layers. The field has a strength 
between 2 and 40 G. Different from the surface magnetic field, it is more ho- 
mogeneous and fills practically the entire space so that the magnetic filling 
factor in the top layer is close to unity. There seems to be no spatial correlation 
between chromospheric flux accumulations and the small-scale field concen- 
trations in the photosphere. Magnetoacoustic waves that form transient fila- 
ments of stronger than average magnetic field are a ubiquitous phenomenon in 
the chromosphere. They form in the compression zone downstream and along 
propagating shock fronts. These magnetic filaments that have a field strength 
rarely exceeding 40 G, rapidly move with the shock fronts and quickly form 
and dissolve with them. Hence, the coupling of waves with the magnetic field 
leads to a continuous agitation of the magnetic field in the chromosphere by 
shock waves. It is not yet clear what the significance and the consequences 
of these perturbations are, especially in view of electro-magnetic dissipation 
processes. 



4 Coupling of radiation with magnetic fields 

The radiative flux in the solar atmosphere couples to the magnetic field 
not only microscopically (changing polarization state) but also macroscopi- 
cally through modification of the gas pressure and density by the magnetic 
field. The magnetic field has an energy density, Cmag, often called the mag- 
netic pressure, like gas pressure has: Pg = ethcrm- The ratio of the two is 
P = ethorm/emag- In locations where /3 < 1, magnetic pressure substitutes 
gas pressure in the transversal direction, which at the same time lowers the 
material density (given thermodynamic equilibrium) and with it the opac- 
ity. Hence, as a rule of thumb, opacity is lower where the magnetic field is 
strong, which leads to a redirection of the radiative flux. This effect was called 
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"flux channeling" by Cannon (19851, because the net radiative flux, which is 
strictly vertical and outwardly directed in a plane parallel stellar atmosphere, 
becomes partially and laterally redirected in the presence of a strong, vertically 
directed magnetic flux concentration. In radiative equilibrium this effect leads 
to a temperature perturbation in the sense that the temperature is slightly 
enhanced in the surface layers of the magnetic flux concentration but lowered 
in the deep layers as is shown in Fig. |8] 




Fig. 8. Temperature contours of a two-dimensional schematic magnetic flux sheet 
in radiative equilibrium with its surroundings. The opacity in the center rectangle is 
one flfth of the opacity in the rest of the domain, causing a redirection (channeling) of 
the radiative flux. In this state of radiative equilibrium the flux channeling causes a 
lowering of the temperature in the bottom part of the "flux sheet" and a temperature 
enhancement in the top part, where radiation starts escaping to the ambient free 
space. Adapted from Steiner ( 1991 1. 



The coupling of the radiation fleld to the magnetic fleld via the channeling 
effect introduces a substantial anisotropy in the radiation fleld at the location 
of a magnetic flux concentration and in its surroundings. This anisotropy gives 
rise to the facular phenomenon. 

Faculae can be explained in terms of the "hot wall model" (Spruit 1976 1. In 
this model the hot wall corresponds to the wall of the depression that is caused 
by a magnetic flux concentration, like the Wilson depression of a sunspot. The 
wall separates the magnetic flux concentration from the surrounding, practi- 
cally fleld-free plasma. However, the term "hot wall" may be misleading as 
it is by far not as hot as the plasma in the same height range in the unper- 
turbed ambient atmosphere. In fact, the rectangular wall of the low-opacity 
region of the radiative equilibrium model shown in Fig. [8] has in most parts a 
lower temperature. In the case of photospheric magnetic flux concentrations 
the temperature of the wall is determined by a delicate balance between ra- 
diative losses and the convective energy supply, where the latter is actually 
reduced in the transverse direction close to the flux concentration because the 
magnetic fleld acts like a solid wall to plasma motion in this direction. Equi- 
librium models and magnetohydrodynamic simulations, however, invariably 
show a bright edge where the wall of the "Wilson depression" , merges with 



the horizontal surface (Carlsson et al. 2004 Keller et al. 2004 Steiner 2005 
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De Pontieu et al.|2006 ) . Seen at an oblique angle (corresponding to an obser- 



vation off disk center) this "hot edge" leads to a characteristic center-to-limb 
variation in contrast. 

It was only recently realized when analyzing high resolution filtergrams of 
faculae from the Swedish Solar Telescope (SST) that their contrast enhance- 
ment extends quite a distance in the limbward direction, typically for half a 
granular size, which is much further than the depth of the depression wall 
would possibly be ( Lites et ar]|2004 Hirzberger fc Wielir||2005 1 . Berger et al. 



( 2007 1 measure an average radial width of faculae of 400 km. This means that 
the contrast enhancement of faculae extends beyond the depression proper 
in the limbward direction. The reason for this behavior is explained with the 
help of Fig. [9] as follows. 




Fig. 9. Photons preferentially escape along 
the line of sight that traverses the magnetic 
flux concentration because of its rarified (less 
opaque) atmosphere. Hence, the radiation field 
lateral to the fiux concentration is asymmetric 
due to radiative flux channeling. 



A material parcel located in the solar atmosphere lateral to a magnetic flux 
concentration "sees" a more transparent atmosphere in the direction toward 
the flux concentration as compared to a direction under equal zenith angle 
but pointing away from it because of its rarefied atmosphere. Consequently, 
from a wide area surrounding the magnetic flux concentration, radiation es- 
capes more easily in the direction towards the flux sheet so that a single flux 
concentration impacts the radiative escape in a cross-sectional area that is 
much wider than the magnetic fleld concentration proper. This means that 
the presence of a magnetic flux concentration introduces an anisotropy in the 
radiation fleld such that when observed close to the limb, a granule limbward 
of (seen across) the flux concentration appears brighter than normal. Hence, a 
facula is not to be identifled with bright plasma that sticks, as the name may 
insinuate, like a torch out of the solar surface. Rather is it the manifestation 
of photospheric granulation, seen across a magnetic flux concentration ~ gran- 
ulation that appears brighter than normal in the form of so called "facular 
granules" . 

Another consequence of the flux channeling effect is an enhancement of 
radiative loss from the solar surface at the location of the small-scale mag- 
netic flux concentration. One could say that magnetic flux concentrations add 
a roughness to the solar surface that increases its effective area and thus in- 
creases the radiative loss from it. Therefore, the coupling between radiation 
and magnetic fleld, and with it faculae, play a key role in the solar irradiance 



variation over a solar cycle and on shorter time scales (Fligge et al. 2000 
Wenzler et ar]|2005| [Foukal et ar]|2006| [Krivova fc Solanki|,200^ 
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5 Conclusion 

Results obtained with the spectropolarimeter of the Solar Optical Telescope 
onboard the Hinode space observatory have greatly extended our knowledge 
of the magnetism of the quiet Sun. It has now become apparent that virtually 
every location on the surface of the Sun harbors magnetic field that can be 
detected in the 630 nm Fel spectral line via the Zeeman effect, exceeding a 
la noise level of 0.6 G and 20 G for the longitudinal and transversal magnetic 
field, respectively ( Lites et al.||2008 l. At the scale of 0.3", this magnetic field 



has a preferential direction, which is parallel to the solar surface, i.e., it is 
anisotropic. It seems now that Hinode has discovered the majority of the 
hitherto "hidden" magnetic field that was known to us only through Hanle 
measurements, and it made it accessible to Zeeman analysis. 

The predominantly horizontal direction of the weak internetwork magnetic 
field is also a robust result of magnetohydrodynamic simulations of the sur- 
face layers of the Sun. It can be explained in terms of the coupling between 
magnetic field and convective plasma motion. Two aspects of this interaction 
are evident. First, the process of flux expulsion, which displaces magnetic field 
from eddy cells. It leads to a concentration of predominantly vertical fields in 
the interganular lanes and predominantly horizontal fields in the stable layers 
above and at the edges immediately adjacent to granules. Second, the turbu- 
lent dynamo, which operates in the top surface layers of the convection zone, 
leads in the adjacent, stably stratified photospheric layer to a multi-scale sys- 
tem of magnetic loops. This loop system naturally leads to a stronger decline 
with height of the vertical than of the horizontal component of the magnetic 
field, hence, to a predominance of the horizontal fields in the height range 
where the spectral lines used for the Hinode observations are formed. It is 
still a matter of future research to find out to which degree the internetwork 
magnetic field is due to the turbulent surface dynamo, the remnants of pre- 
existing magnetic fiux of active regions, and the emergence of magnetic flux 
from the deep convection zone, or due to yet other, additional sources. Also 
the role of the horizontal field in the heating of the chromosphere, e.g., by 
Ohmic dissipation of associated current sheets, needs yet to be clarified. 

The study of the propagation of magnetoacoustic waves in a magnetically 
structured, stratified atmosphere is a relatively new field of research. It has 
applications in magnetoatmospheric seismology, for example, when determin- 
ing the topography of the magnetic canopy by measuring wave travel times in 
the photosphere and chromosphere. The surface defined by equality between 
sound speed and Alfven speed is a zone of mode conversion. It separates the 
magnetically dominated tenuous region of fast moving magnetic modes and 
magnetic shock waves from the more slowly evolving atmosphere beneath it. 
The strong gradients in the Alfven speed of a magnetically structured at- 
mosphere lead to refraction and refiection of the magnetically driven modes. 
Future research in this field is directed at improving diagnostics for magne- 
toatmospheric seismology. Measurement and theory of magnetoatmospheric 
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wave propagation need to be improved and linked for a reliable interpretation 

and exploitation of observations. The role of wave mode conversion and the 
channeling of slow modes in magnetic flux concentrations for the heating of 
the outer atmosphere must yet be quantified. 

Magnetic flux concentrations lead to radiative flux channeling and in- 
creased anisotropy in the radiation field. Faculae can be understood in terms 
of this effect, which is apparent from magnetohydrodynamic simulations in 
two and three spatial dimensions. Future research in this field should include 
a thorough statistical analysis of (synthetic) faculae of three-dimensional sim- 
ulations in order to understand the role of magnetic field stength, size, shape, 
etc., for the center-to-limb behavior of the facular contrast and for the energy 
balance of faculae. The latter is crucial for a better understanding of solar 
radiance variability. 
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